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We utilize accurate experimental data available in the literature to yield bounds on the polar- 
izabilities of the ground and first excited states of atomic Yb. For the 6s 2 So ground state, we 
find the polarizability a to be constrained to 134.4 < a < 144.2 in atomic units, while for the 
6s6p 3 Po excited state we find 280.1 < a < 289.9. The uncertainty in each of these values is 1.0. 
These constraints provide a valuable check for ab initio and semi-empirical methods used to compute 
polarizabilities and other related properties in Yb. 
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I. INTRODUCTION 



Lifetimes, polarizabilities, and long-range interaction 
parameters are intimately connected atomic properties. 
For a particular atomic state, these properties may be 
precisely determined given a complete knowledge of oscil- 
lator strengths, together with associated frequencies, for 
all transitions connecting to that state (we restrict our 
attention to transitions and interactions of the electric 
dipole type throughout). While such a complete char- 
acterization of oscillator strengths eludes experimental 
determination, the atomic properties themselves may be 
measured directly. Knowledge of one or more proper- 
ties (e.g., lifetimes) may be exploited to obtain useful 
information on the other related properties (e.g., polar- 
izabilities or van der Waals coefficients). 

In this paper, we employ accurate experimental data 
from the literature to set upper and lower bounds on 
the polarizabilities of the 6s 2 x 5 and 6s6p 3 Pq states 
of atomic Yb, for which experimental values are lacking. 
Three key parameters we utilize are the 6s6p P° radia- 
tive lifetime, known to better than 0.1% from the work 
of Takasu et al. Q, the 6s 2 %— 6s 2 % van der Waals 
coefficient Cq, known to better than 2% from the work 
of Kitagawa et al. @, and the 6s 2 %— 6s6p 3 P ° differ- 
ential polarizability, known to 0.002% from the work of 
Sherman et al. 3 . We further supplement this data with 
experimental lifetime results compiled in Refs. [1, At 
the la confidence level, our results constrain both po- 
larizabilities to a window of width 12 atomic units. For 
the 6s 2 \Sy polarizability, this is comparable to the best 
theoretical results, while for the 6s6p 3 Pq polarizability, 
this window is half that of the best theoretical results. 

We note the similarity of this work to Ref. [(| , wherein 
the lifetimes of the first two excited states of atomic Cs 
were extracted from the measured ground state van der 
Waals coefficient C§ and ratio of the two lifetimes. An 
essential difference between Ref. @ and the present work 
is that, whereas Ref. @ included additional input from 
sophisticated many-body calculations, our present results 
are based entirely on experimental data, free from all but 
first-principle theoretical input. 



II. PRELIMINARIES 

In this section, we present formulae for the ac polar- 
izability a(co), the atom- wall interaction coefficient C3, 
and the van der Waals coefficient Cq. The expressions 
here assume atoms in the ground state, with atomic units 
being used throughout. 

The ac polarizability a{ui) describes the atomic re- 
sponse to a harmonic electric field oscillating at (angular) 
frequency uj. It may be written as 
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where the summation is over all transitions from the 
ground state, with integration over continuum states be- 
ing implicit. Here /; and a;, represent the oscillator 
strength and frequency associated with the i-th transi- 
tion, respectively (fi,u>i > 0). The static polarizability, 
which gives the response to a static electric field, is ob- 
tained by evaluating the ac polarizability in the zero fre- 
quency limit, a = a(0). The term polarizability, when 
given without specification, will imply the static polariz- 
ability. 

The atom-wall interaction constant C3 describes the 
interaction of an atom with a perfectly conductive sur- 
face. C3 may be concisely expressed in terms of the ac 
polarizability as Q 

1 f°° 
C3 = — / a(iu))du. 
4tt J a 

Note that a(iu) is real; inclusion of the imaginary factor 

i in the argument merely effects a change of sign ► + 

in the denominator of expression ([1]) for the ac polariz- 
ability. With Eq. ([T]), C3 may be written in terms of 
contributions from individual transitions, 

1 

where the integration over lu has been performed analyt- 
ically. 

The van der Waals coefficient Cq, which describes the 
long-range interaction between two atoms, may also be 
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expressed in terms of the ac polarizability Q , 
3 



C fi 



\a(iuj)\ dw. 



Using Eq. (TTJ) and performing the integration over oj an- 
alytically, we arrive at the expression 
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Note that the denominator here prohibits factorization 
of the summations over the indices i and j . 

Another useful property is the summation over all os- 
cillator strengths, J^/*. In the non-rclativistic limit, 
the well-known Thomas-Rcichc-Kuhn (TRK) sum rule 
asserts that this summation is equivalent to the total 
number of atomic electrons Af, i.e., 



Rclativistic effects, however, lead to a departure from 
this simple interpretation; these effects will be discussed 
more below. 

The remainder of this section is dedicated to introduc- 
ing notations and conventions to be used in the following 
sections. We start with the definition an = fi/uif, such 
that an represents a partial contribution to the static po- 
larizability from the i-th transition, 
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It should be noted that these contributions are positive, 
i.e., a,; > for all i. Below, expressions for C3, Cg, and 
Af will be written in terms of on in favor of oscillator 
strengths /j. 

The summations over indices i and j in the above 
equations run over all allowed transitions. We will find 
it advantageous to partition the spectrum into a subset 
comprised of select lower- lying ( "main" ) transitions and 
another subset containing all remaining transitions. We 
use the convention of reserving indices a and b for the 
main transitions and indices r and s for the remaining 
transitions. With this convention, properties a, C3, and 
Af are partitioned into two terms, 



a = a=n + ay, 

a r 

C 3 = ^a,w„+^ 



a r uj r , 



a r ui r 



We will refer to the respective terms on each line as the 
"main" and "tail" terms, i.e., 



„ rt main , ^ tail 

a = a + a , 

£fmain _|_ £ftail 

A/- = A/-main + A/-tail_ 
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The Cq coefficient, on the other hand, decomposes into 
three terms, 
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which we will refer to as the "main," "tail," and "cross" 
terms, respectively, i.e., 

C 6 = C™ ain + C^ ail + C£ ross . (3) 

Finally, we introduce the frequency lu , which we set 
equal to the smallest transition frequency outside of the 
main transitions. 



III. LIMITS ON THE POLARIZABILITY 

Partitioning of the spectrum is motivated by the fact 
that transition frequencies and oscillator strengths for 
several low-lying transitions are known or can be deter- 
mined sufficiently well from experimental data. With 
these transitions comprising the main subset, all "main" 
terms in Eqs. ([2]) and Eq. ([3]) can be determined by di- 
rectly adding the contributions from these transitions. 
Limited experimental information for the remaining tran- 
sitions prohibits such direct evaluation of the "tail" and 
"cross" terms. However, given a Cq coefficient and to- 
tal number of electrons Af, we may set bounds on these 
residual terms. In particular, below we derive upper and 
lower bounds on a tal1 in terms of Cq , Af, and the various 
"main" terms. Constraints on a itself are then simply 
obtained by adding a main to each of these bounds. 



A. Upper bound on a 

To derive an upper bound on a, we begin by consider- 
ing the factor LOiUjj /(uji+ujj), which appears in the expres- 
sion for the Cq coefficient. By taking partial derivatives, 
we find that this factor increases monotonically with re- 
spect to both u>i and uij. With this insight, we establish 
the following inequalities, 
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These inequalities may be used to give lower bounds on 



C*^ ail and Cg ross in terms of a tail , 
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where we have introduced the factor £ according to the 
relation, 
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£ is positive; furthermore, if wo is larger than all main 
transition frequencies, it follows that £ is necessarily less 
than unity as well. 

The bounds on Cg al1 and Cg ross give a corresponding 
bound on Cq itself, 

Q Q 
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Noting that all factors here are necessarily non-negative, 
we may rearrange this inequality to yield an upper bound 
on a tal1 , 
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By simply adding a mam , we arrive at the result 
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a < a mai„ ^ a main + (£ a main)2 + _U 6 >_ _ 

V 3 w 

(4) 

In the limit when no transiti ons are inclu ded within the 
main subset, we find a < y/ (4/3)C6/o;o, with ujq being 
the first allowed transition frequency. 



B. Lower bound on a 

To derive a lower bound on a, we again start with the 
factor ujiUJj/{uJi + ujj), this time utilizing the two inequal- 
ities 



U r UJ s 1 

; < t K + uj s ) , 



(5) 



The first inequality follows from 4uj r uj s < \uj t lo s + 
(uj r — ui s ) 2 = (u! r +u! s ) 2 , while the second is obtained 
by taking the limit uj r — > oo. These inequalities may be 
used to set upper bounds on Cg al1 and Cg ross in terms of 
a taU , 

Cf 11 <^Y1 °r a '\ (<* + w.) = 6Cf u a tail , 
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The equality for Cg ross here only holds if either all or 
none of the transitions are included in the main subset. 
Together these inequalities give an upper bound on Ce 
itself. 
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In absence of direct experimental information for C3 , 
expression ([6]) is of limited utility in its present form. 
However, with inspiration from elementary probabil- 
ity and statistics, we may find a useful limit on Cg 3,1 '. 
Namely, for p r > and J2 r Pr = 1' tnc following rela- 
tions hold, 



where (x n ) = J2 r Pr x r- ^ n the language of probability 
and statistics, this merely states that variance is non- 
negative. By making the associations p r — > a r /a tal1 and 
x r — > io r , we arrive at the inequality, 



C 



tail 



< 



Taking this result together with (J6)), we find 



C 6 < Q nain + -y/Ar-Af maiu (a tail ) ' + 24C ; 
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where we have further used A" tail = J\f- W main . 

Our aim now is to rearrange inequality ([7]) to obtain 
bounds on a . To this end, we note that in the limit 
of equality, expression (JT]) becomes a cubic equation for 
V a tal1 . A general cubic equation ax 3 + bx 2 + cx + d = 
has three solutions for x; for the special case of a, b > 0, 
c = 0, and d < — such as the present case — only one 
solution is a non-negative real number. This solution 
corresponds to a lower bound on V a tai1 . Squaring this 
solution to the cubic equation and adding a mam , we find 
a lower bound on a, 
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where the function is defined for x > by 



2 cosh I — cosh 



1-x 
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The function g{x) increases monotonically with x, ap- 
proaching zero [specifically, g(x) — > 3(x/2) 1 / 3 ] in the 
limit x — > and unity in the limit x — > 00. An equiva- 
lent function may be obtained by simultaneously replac- 
ing the hyperbolic cosine function and its inverse by their 
trigonometric counterparts (i.e., cosh — > cos, cosh -1 — s- 
cos -1 ) Q. In the limit of no transitions being included 

r -1 2/3 

in the main subset, we find a > (4/3)C , 6/vA/' 



IV. EXPERIMENTAL BOUNDS ON THE 6s 2 ^0 
AND QsQp 3 P ° POLARIZABILITIES IN YB 

Inequalities (j4]) and ([8]) represent the principle results 
of the previous section. Here we use these inequalities 
together with experimental data to set constraints on 
the 6s 2 ^Sq ground state polarizability of Yb. Moreover, 
we further constrain the polarizability of the 6s6p 3 Pq 
excited state by employing the result of a recent high- 
accuracy measurement of the 6s 2 — 6s6p 3 Pq dif- 
ferential polarizability. In the following section, these 
constraints are compared with theoretical values of the 
6s 2 ^Sq and 6s6p 3 Pq polarizabilities given the literature. 
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TABLE I. Transition frequencies cj; and contributions to po- 
larizability a, for the ground state of Yb. The lowest six 
transitions from the ground state are displayed, with i en- 
compassing all degenerate magnetic states of the upper level. 
Configurations given with three electrons imply excitation of 
an electron from the otherwise filled 4/ shell. All values are 
in atomic units. 



(6s 2 % ->■ ) 



6s6p 3 Pi 
6s6p 1 Pf 
5d6s 
5d6s 2 I 
6s7p 3 Px° 
5d6s 2 3 Dl 



1 5 1 

2 ' 1i 



0.081978 
0.114219 
0.131482 
0.170473 
0.173934 
0.175065 



} 



2.4±0.„ 
100.4 ±0.1™ 
21.1 ±1.^ 

1.6±0.8 a ' b 




A. Bounds on the 6s 2 1 S polarizability 

In Table [IJ we compile data for the lowest transitions 
from the ground state in Yb. Transition frequencies 
are taken from Ref. Q, while partial contributions ct!j 
are inferred from experimental lifetime data given in 
Refs. [H, 0, H[. We ascribe uncertainties to a.% based 
on uncertainties quoted in the relevant references. The 
6s 2 x Sq — > 6s6p X P^ transition gives the dominant contri- 
bution to the polarizability (~70%). With an accurate 
value for the 6s6p 1 P? lifetime known from photoassoci- 
ation spectroscopy [if, this dominant contribution is de- 
termined to within 0.1%. The 6s 2 % -> 5d6s 2 (f , |)° 
transition, involving an excitation from the 4/ electronic 
shell, gives the next largest contribution to the polariz- 
ability (~ 15%). We have determined this contribution 
by taking the weighted mean of five lifetime results com- 
piled in Ref. Q • Rather than using the conventional un- 
certainty associated with the weighted mean, 0.7 in this 
case, we adopt a more conservative value of 1.2, which is 
commensurate with the smallest uncertainty of the five 
lifetime results. 

The last three transitions in Table U have similar fre- 
quencies. The proximity of these three states to one an- 
other is cause for some special consideration. For ex- 
ample, based on configuration and spin multiplicity, one 
may suspect the first of these states to decay almost ex- 
clusively to the ground 6s 2 ^Sq state, with the second 
decaying almost exclusively to the 5d6s 3 D 1 and 6s7s 3 S 1 
states. The validity of configuration and spin quan- 
tum numbers, however, may breakdown due to strong 
Coulomb mixing between these three states. While the 
near-degeneracy brings this additional complication, it 
also allows for a counterbalancing simplification by per- 
mitting us to consider the three states simultaneously. 
Based on experimental lifetimes from Refs. [HQ, we take 
the cumulative contribution to the polarizability from 
these states to be 1.6 ± 0.8; the 50% uncertainty here is 



a conservative estimate of the potential effects of mixing 
between these states. We take a single transition fre- 
quency LJi to be representative of these transitions; our 
results below are insensitive to the particular choice of 
uii when taken within the window of the three transition 
frequencies. 

The next transition from the ground state, beyond 
those shown in Table [I] is the 6s 2 ^Sq — ► 6s7p X P° transi- 
tion at uii = 0.184823. Experimental lifetime results from 
Ref. [3] indicate that this transition could provide a con- 
tribution as large as ~ 9 to the polarizability. This would 
be the case if the 6s7p X P^ state decayed exclusively to 
the ground state (an assumption which appears to be 
made in Ref. [HI). However, electric dipole coupling to 
the 5d6s 1 D 2 and 6s7s S states may also be strong. 
Lacking essential information on the branching ratios of 
the 6s7p P° decay, we opt not to designate a value 
for this transition. Consequently, we take the transitions 
in Table U to comprise our main subset, and assign to ojq 
the value of the 6s 2 ^ — > 6s7p X P° transition frequency. 

With eti given for all transitions in the main subset 
(Table U, the factors a main , £a main , Cf ain , C^ ain , and 
Af main may be determined directly. To obtain bounds on 
a from inequalities ((4]) and ((U) , we further require the van 
der Waals coefficient Cq and total number of electrons Af. 
Kitagawa et al. have probed the long-range interac- 
tion between ground state Yb atoms via photoassociation 
spectroscopy; in their work, the authors determined the 
van der Waal's coefficient to be 



C R = 1932 ± 30. 



(9) 



Taking contributions from Table [J together with this 
Cq coefficient and Af — 70, we obtain the following 
bounds on the ground state polarizability of Yb, 

134.4 ± 1.0 < a < 144.2 ± 1.0 . 

The uncertainties given here were obtained by performing 
Monte Carlo calculations of the bounds [r.h.s. of inequali- 
ties (|2]) and ijHJ)] starting with normally distributed values 
of (Xi from Table U and Cr from Eq. ((9]). We note that the 
uncertainty in each bound, 1.0, is below the uncertainty 
of a mam itself (sec Table []}; this is due to covariance be- 
tween the "main" terms which appear in inequalities (j4]) 
and ([!]). 

We recall that the number of electrons Af, which we 
used here to obtain a lower bound on a, appeared in our 
expressions by invoking the TRK sum rule. The TRK 
sum rule, however, is only strictly valid in the nonrcla- 
tivistic limit, and appreciable deviations might be sus- 
pected for a heavy atom such as Yb. The lowest or- 
der relativistic corrections scale as (Z/c) 2 , with Z be- 
ing the nuclear charge {Z = Af for a neutral atom) and 
c 137 being the speed of light. We find that making 
the substitution Af -> Af [l ± {Z/c) 2 ] = 70 ± 18 in ex- 
pression ((HJ only shifts the lower bound by =F0-3, this 
being well within our uncertainty. Moreover, dedicated 
calculations have found the relativistic corrections to the 
TRK sum rule to be only «1% for Yb [ll|, suggesting 
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that the above substitution grossly overestimates these 
effects. Therefore, we conclude that it is acceptable to 
neglect relativistic corrections to the TRK sum rule in 
the present case. 



B. Bounds on the 6s6p Po polarizability 

The atomic structure of Yb is well-suited for a fre- 
quency standard based on neutral atoms confined in an 
optical lattice trap [HI, [l3[ . Over the course of the last 
few years, Yb lattice clocks have demonstrated perfor- 
mance on par with the best neutral atom and single-ion 
frequency standards, while still holding potential for fur- 
ther improvement [l5j]. One of the largest system- 
atic shifts to the clock frequency is due to thermal ra- 
diation from the room temperature environment imping- 
ing upon the atomic sample [ltjj . In an effort to reduce 
the uncertainty associated with this thermal shift, Sher- 
man et al. Q recently measured the static Stark shift to 
the 6s 2 1 S 6s6p 3 Pq optical clock frequency to high 
accuracy. Their results were reported in terms of the 
differential polarizability, 



a (6s6p 3 P ° 



(6s 2 % 



145.726 ±0.003 



We use this result, together with our results above for 
the ground state polarizability, to set constraints on the 
polarizability of the 6s6p 3 Pq excited state, 

280.1 ± 1.0 < a < 289.9 ± 1.0 . 

The excited state polarizability is about twice as large 
as the ground state polarizability; as a consequence, the 
constraint on the excited state polarizability is fraction- 
ally about a factor of two more narrow than the con- 
straint on the ground state polarizability. 



V. COMPARISON WITH THEORETICAL 
VALUES FROM THE LITERATURE 

In Figure [TJ we display results of ab initio and semi- 
empirical calculations taken from the literature for the 
6s 2 and 6s6p 3 Pq polarizabilitics. Along with these 
values, we further display the present constraints derived 
from experimental data. 

For the ground state, we find that several of the the- 
oretical values lie within our constraints. In particular, 
with the exception of Ref. [l6| , all results with explicit er- 
ror bars arc found to agree very well with our constraints. 
Dzuba and Derevianko Q identified an oversight in the 
semi- empirical method of Ref. (l6j . which accounts for 
the discrepancy with this result. At the la confidence 
level, we see that our present results and the most accu- 
rate theoretical results of Refs. @, H3] constrain the po- 
larizability to a similar window, having a width of about 
12. 

For the excited state, we find that the theoretical re- 
sults with explicit error bars agree reasonably well with 
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FIG. 1. (color online) Theoretical values from the literature 
for the polarizabilities of the 6s 2 1 Sq (top panel) and 6s6p 3 Pq 
(bottom panel) states of Yb, together with the present con- 
straints derived from experimental data. Theoretical values 
are denoted with circle markers, along with error bars where 
available. The vertical bands indicate the present constraints, 
with the left band extending la below the lower bound and 
right band extending la above the upper bound. 
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our constraints. Specifically, values from Refs. [l6l . l2Cil . l26j 
are seen to lie about la below our lower bound, while the 
value from Ref. Q lies about la above our upper bound. 
The only other result, having no explicit error bar, also 
lies above our upper bound. At the la confidence level, 
we see that our results constrain the polarizability to a 
tighter window — about a factor of two narrower — than 
the best theoretical results. 

It is not the goal of the present paper to provide criti- 
cism of the theoretical methods used in Refs. 0, [Mill- 
ie . Nevertheless, as a general comment, it appears that 
our constraints give a certain degree of validation for sev- 
eral of these methods. We anticipate that our results will 
serve as a valuable check for ab initio and semi-empirical 
methods used in future works to calculate polarizabili- 
ties, as well as other related properties of Yb (e.g., C 3 
and Cq coefficients, magic wavelengths, dipole matrix el- 
ements). Fig. [1] illustrates that our present constraints 
are relevant for this purpose. 



VI. COMMENTS ON THE SPECTRAL 
DISTRIBUTION OF a tail 

The gap between our upper and lower bounds demon- 
strates our ignorance in the spectral distribution of a tal1 
across the frequencies u> r > u>$. For example, the equal- 
ity of @, which gives the upper bound on the polar- 
izability, is only satisfied if the entirety of a tal1 is ac- 
cumulated from transitions with frequency precisely lvq 
and no larger. With an appropriate model for the dis- 
tribution of a tal1 (or, cquivalently, distribution of oscil- 
lator strengths), we could arrive at a central value for 
the polarizability lying somewhere between our present 
constraints. We have chosen not pursue such a program, 
stressing that our present results are derived from exper- 
imental data available in the literature and are indepen- 
dent of such theoretical modeling. 

Unlike upper bound (£[]), the lower bound (J5J does not 
correspond to a precise distribution of a tal1 . From nu- 
merical calculations, we found that, with Cg, J\f, and 
"main" terms fixed, a tal1 is minimized by a distribution 
peaked sharply about a frequency lo t « 2.5. This real- 
ization motivated the choice of inequalities in (|5|). The 
first is near equality when \u) r — lu s \ -C (co r + oj s ), while 
the second is near equality for tu a <C oj r . Therefore, with 
these inequalities our analytical lower bound (0) has a 
close correspondence to this minimal distribution. While 
it is physically unlikely that a tal1 is accumulated entirely 
at frequency u> r w 2.5, we again emphasize that we do 
not speculate on the actual spectral distribution of a tal1 . 



VII. EXTENSION TO NON-STATIC 
POLARIZABILITIES 

While the primary focus of this work is to obtain con- 
straints on the static polarizability a = a(0), the results 



above can further be extended to yield constraints on 
the ac polarizability a(oj) evaluated at other frequencies 
u 7^ 0. To illustrate this, we consider the ac polarizability 
evaluated at the "magic" lattice frequency u* , which bal- 
ances the ac polarizabilities of the 6s 2 ^Sq and 6s6p 3 Pg 
clock states [12|. Dedicated measurements have deter- 
mined the magic frequency to be w* = 0.06000 fl5l. [29j. 
Defining 8 = a(u)*), it follows from Eq. (jTJ that contri- 
butions to 8 satisfy 

3 a ^ 

Ujf — LJ* Z 

The "main" term /3 mam = J2 a Pa c an be computed di- 
rectly from the values given in Table U while the "tail" 
term /3 tal1 = ^2 r Pr is necessarily limited to the range 

„ tail ^ otail ^ „ tail ( ^0 \ ( ~\ c\\ 

a </3 - a Uf^jJ' (10) 

where the factor in parenthesis evaluates to 1.12. Using 
inequalities (fTOj) together with inequalities (|4|) and © for 
bounds on a tal1 , we obtain the following constraints on 

^ #main _|_ ^tail 

181.1 ±1.3 </3< 193.2 ±1.1 . 

These constraints have been derived specifically for 
the 6s 2 Sq ground state; however, we note that 
they are equally applicable for the 6s6p 3 Pq excited 
state, following from the definition of magic frequency, 
i.e. 8 (6s6p 3 P °) = B (6s 2 %). Knowledge of 8 = a{to*) 
is useful, as it can be used to directly relate lattice inten- 
sity to trap depth in the optical lattice clock j29j . 

VIII. CONCLUSION 

We have utilized experimental data available in the 
literature to yield bounds on the polarizabilities of the 
6s 2 ^ and 6s6p 3 P ° states of Yb. Key experimental 
parameters employed were the 6s6p P° radiative life- 
time from Takasu et al. [l[, the 6s 2 ^Sq — 6s 2 ^Sq van der 
Waals coefficient Cq from Kitagawa et al. Q, and the 
6s 2 ^Sq — 6s6p 3 Pq differential polarizability from Sher- 
man et al. Q. For the 6s 2 state, our results constrain 
the polarizability to a window comparable to that of the 
best theoretical results (at the la confidence level), while 
for the 6s6p 3 Pq state, our results constrain the polariz- 
ability to a window half the width of the best theoreti- 
cal results. We anticipate that our results will serve as 
a valuable check for ab initio and semi-empirical meth- 
ods aimed at calculating polarizabilities and other related 
properties in Yb. 
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